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We study theoretically the optics in undoped direct gap semiconductors which are strongly driven 
in the THz regime. We calculate the optical sideband generation due to nonlinear mixing of the THz 
field and the near infrared probe. Starting with an inversion symmetric microscopic Hamiltonian 
Ch , we include the THz field nonperturbatively using non-equilibrium Green function techniques. We 

find that a self induced relativistic spin- THz field coupling locally breaks the inversion symmetry, 
resulting in the formation of odd sidebands which otherwise are absent. 
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rH ■ I. INTRODUCTION 

^ . . 

O \ Optical properties of semiconductors are sensitive to external conditions such as strong static electric fields being 

^ applied to the sample. This was predicted more than 40 years ago by Franz and Keldysh ^ and leads to the Franz- 

^ ; Keldysh effect. The effect manifests itself in finite absorption within the gap, near the band edge, and modulation 

d . of the above gap absorption spectrum. The effect can be understood in terms of tunneling assisted absorption. This 

^ ■ work was generalized to oscillating external fields, E{t) = i?cos(f2t), by Yacoby |^ 30 years ago. Yacoby concluded 

I , that similar effects to the Franz-Keldysh effect manifest provided that the ponderomotive energy, 
X3 ■ 

Q : Ef^nu^f = — — , (1) 
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magnitude is similar to the photon energy, 7if2, of the external field, i.e. 



Here l/rrir = l/rric + l/rriy is the reduced mass of the effective valence and conduction band masses rriy and rric 
respectively. Then the field induced effects extend a few around gap and may be viewed as being due to a 
■ combination of quantum mechanical tunneling and multi photon processes. Experimentally it is exceedingly difficult 
! to fulfill the condition (|^) such that the effect extends over a range which can be resolved experimentally. The advent 
^ ' of the free electron laser as a source of intense coherent radiation in the THz regime |^ has made it possible to enter the 
regime 7 '--^ 1, which now has come to be known as the dynamical Franz-Keldysh (DFK) regime. Works reporting on 
^ optical experiments in the DFK reg ime include Rcfs. [ This has also lead to renewed theoretical interest of THz 

electro optics, see e.g. Refs. [ |9 -13 ]. It has been shown that the physical signature of the DFK effect is enhanced with 



reduced dimensionality |10[| , which has lead to experimental investigation of the two dimensional analogue of Yacoby's 
prediction in quantum wells at low temperatures Q. Of particular relevance to the present work is the observation 
^ , of optical sidebands in the transmission spectrum at frequencies LOp ± 2nn, where Up is the probe frequency and n 
IL* ' is an integer. Thus the sidebands only appear at even multiples of the driving frequency ft in quantum wells. This 
. !^ I observation can be understood as being a direct consequence of the underlying inversion symmetry of the system 
10, e.g. the dispersion has the property e„(fc) = e„(— fc), n is the band index. We remark from the onset that it is 
, assumed that the THz field does not induce interband coupling ||l^ . Absence of odd sidebands has also been observed 
- - 1 in the presence of strong quantizing magnetic fields j^. In Ref. [ it has now been experimentally demonstrated 
that odd sidebands appear if the inversion symmetry is broken by driving the system such that the field oscillates 
in the growth direction of an asymmetric quantum well. These observations support the assumption that the THz 
field only induces intraband dynamics. In the light of these observations and the theoretical understanding thereof it 
therefore comes as a curious surprise that odd sideband formation is observed of resonance, e.g. away from the gap, 
in apparently inversion symmetric bulk samples [ p^ . In the present work we theoretically demonstrate that inversion 
symmetry breaking results from spin-dynamic field coupling leading to the formation of odd sidebands. We refer to 
the well known relativistic effect that a charged particle traveling with velocity v in an electric field E experiences 
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an effective magnetic field proportional to v x E which couples to the particles spin degree of freedom via a Zeeman 
term, see e.g. Ref. [ [^]. In the present case the electric field is the strong THz field driving the system. 

The article is organized as follows; in Sec. |l| we define the theoretical model and calculate the sideband intensities 
neglecting the spin-THz field coupling and show how only even sidebands appear, in Sec. ^ we calculate effect of 
the spin-THz field coupling and finally in Sec. ^ we discuss the physical consequences of our results and present 
numerical calculations for sideband generation in GaAs and InAs. 



II. THE MODEL 



In this section we define the model which we shall study. Our approach is to apply non-equilibrium Green function 
techniques to include the THz field nonperturbatively. This can be done analytically and the resulting Green functions 
form the starting point for our subsequent calculations, which is to account for the relativistic term perturbatively 
in the Born approximation. The resulting Green functions are then applied to determine the intcrband susceptibility 
which describes the linear response of the system to a weak optical near infrared probe field. 

We consider a direct gap undoped bulk semiconductor subject to an intense linearly polarized THz field, ETHz{t)- 
We shall assume that effects due to the finite wave-vector are negligible, that is we assume that the THz field is 
uniform. Thus the system considered remains translationally invariant and the wave vector is still a good quantum 
number. We shall model the system with a four band model which is a generalization of the model introduced in Ref. 
[ jist]. A conduction band and a valence band, each carrying one of two spin states. We start our analysis with the 
second quantized Hamiltonian 

i?°= E (3) 

Here a is the spin state index, a the band index and k is the wave vector of the state. The THz field is introduced 
nonperturbatively via the vector potential A — —Esvn.{U,t)/^. We shall assume that the THz field is oriented in the 
z-direction. The probe has wave vector q and is oriented such that it forms an angle 9 to the i-direction, see Fig. 
For simplicity, we shall assume that the bands are parabolic, even away from the gap. We assume that 

= ^ + (4) 

for the conduction band. Here TOc is effective mass for the conduction band and eg is energy gap of the semiconductor. 
For the valence band we describe the dispersion as 

where the m„ is the effective mass of valence band electrons. Note that eaa{p) — £aa{—p)i so the initial Hamiltonian 
is inversion symmetric. Furthermore, we shall assume that the interband dipole matrix elements dcva^a^ are locally 
independent of the wave vector, but that they depend on the general underlying symmetries of the Brillouin zone 
being probed. For instance, near a point of high symmetry like the F-point we take dcvu^a^ oc 5a^a^ , reflecting the 
near gap selection rules. Away from the F-point these selection rules are not in force, and we shall assume that all 
the matrix elements are finite but constant. Below we shall thus consider two distinct regimes, the near gap regime 
characterized by the selection rules and the far regime, where the selection rules are not in force. 

The relevant quantity to study in order to glean the optical properties of the system, is the interband susceptibility. 
The non-equilibrium causal interband susceptibility on the Keldysh contour is given by 

i f dk 

Xcvf3q,t,t') = J j^--:^gca,„,^ik + q,t,t')gy„,^„^{k,t',t), (6) 
where (3 — {ac'^ya'^a'y) and the contour ordered Keldysh Green functions are defined by 

gc.aa'{k,t,t') - -*(T,[c_,Kt)ct,^,-(t')]). (7) 

For an overview on non-equilibrium Green functions see p^ . The physically relevant susceptibility is the real time 
retarded susceptibility. We analytically continue the causal susceptibility using the Langreth rules and find that 
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i f dk \ " -> -> -> 1 

XcvffiQAt') = -j-J j^[9c.^a'Sk + qAt')g:^-^^k,t\t)+g:,^,,ik + q^^^^^ (8) 

In the quasiparticle picture one describes the lesser function, g"^, as a distribution function / times the spectral 
function via gf^ =ifa. The spectral function is given hj a = i{g'^ — g"') in a non-equilibrium system. In equilibrium 
this gcncrahzation becomes the familiar a^q = 2Img'^. 

In an undoped semiconductor all the valence band is occupied, hence the distribution function is one, but no states 
in the conduction band are occupied and the corresponding distribution function is zero. We thus put g^ = —{gl — g") 
and g^ =0. The inter band susceptibility becomes 

i f dk 

XlvM^.t') = -^J (2^5:-c<(^ + q,t,t')g:^,^^ik,t',t). (9) 
Probing the system with a weak beam Ep{q,Wpt) E'pe""^"*, the linearly induced interband polarization is 



Driving the system with THz frequency fl the induced polarization will be of the form 

Pe.(t) = e2£;p^r,„(u;p)e-'(--+"")*, (11) 

n 

which spectrally is a comb of oscillating dipoles giving rise to sidebands irradiating with intensity proportional to 

In{^p) = \Vn{^p)\'^ and frequency uj.p + nO.. The absorption of the probe is proportional to lm.rio{u>p). 
Neglecting relativistic effects, the Green functions obeys the Dyson equation 

{ihdt - e„,{nk + eA{t))}g°jJ^,{k,t,t') = h6„„,6{t - t'), (12) 

with the appropriate boundary conditions. These Green functions which include the THz field to all orders form the 
starting point of our subsequent calculations. Thus our unperturbed Green functions take the THz field into account 
from the onset. The equations are readily integrated with the results 

=±?i^<.<.'^(±tTOexp{ -iy^'ye„<.(?ifc + eA(s))}. (13) 

Using these Green functions in order to evaluate the interband susceptibility leads to a Gaussian integral in fc-space 
which we evaluate with the result 

Here t = t — t', T = {t + t')/2, M = my + rric-, and we have defined 

w/ ^ 0{x) /,. sin^(a;/2)\ ^ ( , . ,sin^(a;/2)^\ 

/"(^) = ^sin'(a;/2)exp(^4i7 LA^j ^(^^(sina; - 4 ^))- (15) 

Other components of the susceptibility are zero. We thus obtain that 

r?2.(c^p) = 24.TT^c.TT (^^^3/2 \^ dT!l(SiT)e^^-\ (16) 

where we have used that dcv\\ = dcvH and we have defined ujqn = + nfi — hq^ j (2M) — w/ — LOg. The remaining 
integral is readily evaluated numerically. 
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A. PHYSICAL IMPLICATIONS 



Here several physical results eraerge. Sidebands do only appear at frequencies Wp ± 2ni7, no sidebands involving 
an odd number of THz photons appear. This is a direct consequence of the inversion symmetry of the Hamiltonian. 
The response is independent on the relative orientation of the probe and the THz field, which is reflected by that the 
result only depends on the magnitude of q. In Fig. || we illustrate the sideband intensities for the first 4 sidebands 
as a function of the probe frequency. These results are valid in the near zone, and are in qualitative agreement with 
the experimental findings reported in [Q for quantum wells. Both the probe and the irradiating sidebands are close 
to real states, hence the asymmetry around ujg. The relation In{^) — I—ni^ + nft) holds. The probe or the sideband 
is near the main spectral feature, the band edge in the present case. In the quantum well the main feature is the Is 
exciton resonance 0. 



III. THE RELATIVISTIC SPIN-THZ FIELD COUPLING 



In this section we shall take into account the relativistic effects. The coupling to the spins comes physically about 
because a moving electron in the presence of an electric field experiences a local magnetic field, proportional to v x E, 
which couples to its spin via a Zeeman interaction, see e.g. |1^. This is the effect which gives rise to spin orbit 
coupling in atoms and solids. In the situation considered here the strong THz field provides the electric field giving 
rise to the effect. We describe the contribution to the Hamiltonian as 



U,{k,t) ^ a^ea ■ {k X E{t)), (17) 

here ai, i — {c, v}, is an effective coupling constant which depends on the material, t? is a vector of the Pauli matrices 
and E{t) is the intense THz field. With E oriented in the z direction ( p7| ) becomes 

U,ik,t) = a,e ^^j^^ ) E^. (18) 

Physically the same kind of term has been studied by Raman spectroscopy in asymmetric GaAs quantum wells [ pO| . 
Then the external field is static and comes about due to the asymmetry of the quantum well confining potential. The 
term breaks time reversal symmetry and thus infiuences weak localization leading to weak anti- localization pl| , p2| . 
The important thing for the present study is that Uc has been determined both experimentally |^ and theoretically 
^l| for GaAs leading to ac ~ SA^, and for InAs etc ~ IIOA^. From the point of view of an effective mass picture we 
estimate the coefficient for the valence band by — acimc/m^Y' . These parameters are such that they lead to small 
corrections compared to the remainder of the Hamiltonian. We thus calculate the correction to the Green functions 
to lowest order in the Born approximation, with the result 



/oc 7 
'^gT'''{k,t,s)U,{ls)i^'\k,s,t') (19) 

_ ( 1 ±k{t,t'){ky + ik^) \ Or/a,r , .1^ 

-\±k{t,t'){ky~lk,) 1 ^^"^ 

where we have introduced the quantity 

h{t, t') - cos (r!^) sm{n{t - t')). (21) 

Note that li{t,t') = —li(t',t). Again in this case it is straightforward to perform the k integration in order to obtain 
the interband susceptibility via Eq. (H). The integrals involved are Gaussian times polynomials in k. We remark 
that higher order contributions preserve this structure as well and are thus readily determined. Now the relative 
orientation of the probe with respect to the driving field becomes important. We assume that q is perpendicular to 
the X direction. We find that XcvM r, T) = x°M<l^ r, T) for (3 e {TTTT, Titi, itit, iiii}, 

(2^)3/2 nn ^ 

(f]T)^e^"'^/V^(^^T)e-'2"f^^, (22) 



X cos 
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for/3e {iiit,TiTT,TTTi,iTii}, 

X™M9>^,Jj- (2^)3/2 ^ me ^ mj' 

n 

for /3 e {iTTT) iitii Ttiti Tiii}- Other components are of second order in a and are negligible. We can now find the 
sideband intensities from 

gix(l-2«)/4 ^^3/2 roo 

r,2n{^p) = 2(4.tt4.tT + dcvtldcvu) ^^^-^3/2 'y^ dT/,'^(OT)e'"-^ (24) 

and 

r?,„+iK)^44.TT4.Ti (2,r)3/2 ' dr(/i(f^r) - */i+,(f7r))e--^ (25) 

This is the main result of the present work. 



IV. RESULTS AND DISCUSSIONS 



When dcv^i dcvw, as is the case near the F-point, near resonance, the relativistic effects are suppressed and 
we obtain the results of the previous section. The relativistic effects, i.e. contributions to the sideband intensity 
proportional to ai, manifest themselves in the far zone, away from resonance. The relativistic term leads to the 
formation of odd sidebands. Thus the inversion symmetry has been broken and relativistic self induced local symmetry 
breaking has taken place. The inversion symmetry breaking is local because it is only detectable with finite probe 
wave vector q. We say that the symmetry breaking is global if it leads to odd sidebands even when q —f 0. An example 
of a global symmetry breaking is for instance when a linear term is present in the dispersion. The nature of the even 
( p^ ) and the odd ( p5| ) sidebands is fundamentally different. Let 9 be the angle between the probe beam, and the 
direction the electric field is polarized in, z, see Fig. |l|. The even sidebands do not depend on 0. They can thus be 
regarded as a result of pure amplitude modulations of the fundamental probe. For the odd sidebands however, we 
note that /2n+i oc sin^ 9. The odd sidebands therefore vanish when the probe is aligned parallel with the electric field, 
but is at a maximum when they are aligned perpendicular to each other. 

Taking dcvn — dc-uTT t)e constant, we have evaluated the sideband intensities according to ( p^ and (|2^) using the 
material parameters for the heavy hole band and conduction band for GaAs and InAs. We have taken hil = 30meV 
for the calculations. In Fig. ^ we show /„, n = 1,2, 3, 4, for THz intensities corresponding to 7 = 0.1, 0.5, 1.0, 2.0 as 
a function of the probe frequency for GaAs. Near resonance I2 is dominating. Moving away from the resonance Ii 
decays much slower than the other sidebands and eventually becomes dominating. Notice also that I4 > I3 within the 
region shown. The even sidebands are more sensitive to the THz intensity than the odd sidebands. In Fig. ^ we show 
the results for InAs. Qualitatively the behavior is identical except that for the higher THz intensities we note that 
I4 < I3 < I2 < Ii away from resonance. Quantitatively the odd sidebands are about a factor 100 stronger than for 
GaAs. In Fig. | we show the same sideband intensities as a function of the THz intensity keeping the probe frequency 
fixed, w/rj = 2.5,5.0,7.5,10.0, for GaAs. Note how I2 and I3 seem to rise in a similar manner. It is noteworthy as 
well how feature less Ii is, in fact Ii is mostly linear 7 within the regimes we have investigated. In Fig. ^ we show 
the same results for InAs. 

We have considered other possible sources of inversion symmetry breaking. Spin-orbit splitting due to the crystal 
field leads to a cubic contribution in the effective Hamiltonian for zinc-blend structures [ p3| , which breaks the inversion 
symmetry. We find however that the contribution to the sidebands due to this term is three orders of magnitude less 
than the contribution due to (|7|) considered above for GaAs and InAs. Another possible source is the band bending 
at the edge of the sample due to the pinning of the Fermi level associated with charge accumulation of residual charge 
carriers at the surface of the sample. Thus the relevant inversion symmetry may be broken if the THz field has an 
electric field component perpendicular to the surface of the sample. We have assumed that this is not the case. 

In summary, we have studied the nonlinear generation of optical sidebands in strongly THz driven undoped direct 
gap semiconductors. We have found that relativistic spin-THz field coupling leads to local breaking of inversion 
symmetry which results in the formation of, otherwise suppressed, odd sidebands in the transmitted wave of a weak 
near infrared interband probe. We have shown that the even sidebands are independent on the relative orientation of 
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FIG. 1. Schematic illustration of the experimental setup, showing also the relative orientation of the probe and the linearly 
polarized THz electric field. 



the probe and the hnearly polarized THz field, while the odd sidebands depend strongly on the relative orientation of 
the probe to the THz field. The even sidebands dominate the transmission spectrum for probe frequencies near the 
gap, but moving the probe away from the gap odd sidebands will eventually dominate. We find that the relativistic 
effects are orders of magnitude stronger in InAs than in GaAs. 
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FIG. 2. Sideband intensities as a function of the prob(> frcxjueney using material parameters for GaAs, 7 = 1 and 
TiQ, = 30meV. Only even sidebands appear due to the inversion symmetry of the Hamiltonian. Note that /n(t^) = I-ni'jJ + nQ,). 




FIG. 3. Sideband intensities In, n € {1,2,3,4}, as a function of the probe frequency using material parameters for GaAs 
including the relativistic effects. Values of 7 are 0.1, 0.5, 1.0 and 2.0 with %Q, = 30meV. 
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FIG. 4. Sideband intensities In, n G {1,2,3,4}, as a function of the probe frequency using material parameters for InAs 
including the relativistic effects. Values of 7 are 0.1, 0.5, 1.0 and 2.0 with = 30meV. 




FIG. 5. Sideband intensities In, n € {1,2,3,4}, as a function of 7 using material parameters for GaAs including the 
relativistic effects. Values of w — Wg are 2.5f2, 5.00, 7.50 and lO.Ofi with HQ = 30meV. 
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FIG. 6. Sideband intensities /„, n G {1,2,3,4}, as a function of 7 using material parameters for InAs including the 



relativistic effects. Values oi u) — u)g are 2.5f2, S.Ofi, 7.5fi and lO.Ofi with = 30meV. 
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